Gamma-rays from ultracompact minihalos: potential constraints on the primordial 

curvature perturbation 
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Ultracompact minihalos (UCMHs) are dense dark matter structures which can form from large 
density perturbations shortly after matter-radiation equality. If dark matter is in the form of 
Weakly Interacting Massive Particles (WIMPs), then UCMHs may be detected via their gamma-ray 
emission. We investigate how the Fermi satelhte could constrain the abundance of UCMHs and place 
limits on the power spectrum of the primordial curvature perturbation. Detection by Fermi would 
put a lower limit on the UCMH halo fraction. The smallest detectable halo fraction, /ucmh ^ 10~^, 
is for MucMH ~ 10^ Mq. If gamma-ray emission from UCMHs is not detected, an upper limit can 
be placed on the halo fraction. The bound is tightest, /ucmh ^ 10~^, for Mucmh ~ IO^Mq. 
The resulting upper limit on the power spectrum of the primordial curvature perturbation in the 
event of non-detection is in the range 

^TC ^ 10"'^-^ - lO"'^ on scales fc ~ 10^ - 10^ Mpc'^ This 
is substantially tighter than the existing constraints from primordial black hole formation on these 
scales, however it assumes that dark matter is in the form of WIMPs and UCMHs are not disrupted 
during the formation of the Milky Way halo. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The power spectrum of the primordial curvature per- 
turbation on large scales has been accurately measured 
using cosmological observations fll| . These measurements 
can be used to constrain models of inflation (see e.g. 
Ref. [2]). Cosmological observations only probe a very 
narrow range of scales, however avoiding the overproduc- 
tion of Primordial Black Holes (PBHs) constrains the pri- 
mordial perturbations over a wide range of smaller scales. 
If the density perturbation at horizon entry in a given 
region exceeds a threshold value, 5^ ^ 0.3, then grav- 
ity overcomes pressure forces and the region collapses to 
form a PBH with mass of order the horizon mass 
The abundance of PBHs formed is constrained by the 
consequences of their evaporation and their gravitational 
effects (for recent coinpilations and updates of the con- 
straints see Refs. [1, [Ij). These abundance constraints 
can be translated into constraints on the power spec- 
trum of the primordial curvature perturbation of order 
Pn < 10-1 - 10-2 g. 

Ricotti & Gould [1] have recently proposed that 
slightly smaller perturbations, in the range lO^"^ to Sc, 
can collapse before z ~ 1000 and seed the formation of 
ultracompact minihalos (UCMHs). Due to their early 
formation, the central regions of UCMHs would have a 
high dark matter (DM) density. If DM is in the form of 
Weakly Interacting Massive Particles (WIMPs), WIMP 
annihilation within UCMHs may lead to an observable 
gamma-ray signal ^, 

Scott & Sivertsson [3| have investigated gamma-ray 
emission from UCMHs formed from perturbations which 
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enter the horizon at three different epochs in the early 
Universe: e+e" annihilation, and the QCD and elec- 
troweak (EW) phase transitions. They find that an 
UCMH corresponding to the e+e" annihilation epoch, 
which has present day mass Mucmh (z = 0) ~ lO^M©, 
could be detected by the Fermi satellite or current Air 
Cherenkov telescopes (ACTs), at a distance of 100 pc. 
If 1% of the DM is in the form of UCMHs with this 
mass there would be ~ 3 UCMHs within 100 pc of the 
Earth [7j. UCMHs formed at earlier epochs would be 
lighter, and hence more challenging to detect. 

We extend this work and examine the constraints 
which would be placed on the primordial curvature per- 
turbation power spectrum by the possible detection of 
UCMHs. It has been shown that there are single field 
models of inflation which are compatible with cosmolog- 
ical observations and where the perturbation amplitude 
on small-scales is large enough to produce a significant 
density of PBHs [s", ¥| (see also references therein). It is 
therefore possible that UCMHs may form from perturba- 
tions generated by single field slow roll inflation. Phase 
transitions 0] or features in the inflationary potential 
TO] could also lead to enhanced perturbations on small 
scales. In tliis paper we do not fix the UCMH mass or 
abundance. Instead we calculate the constraints on the 
UCMH halo fraction which would arise from the detec- 
tion (or non-detection) of gamma-rays from UCMHs by 
Fermi as a function of UCMH mass. We then trans- 
late the UCMH abundance constraints into constraints 
on the power spectrum of the primordial curvature per- 
turbation, as a function of scale. In Sec.|TT]we summarize 
the calculation of the properties of the UCMHs and the 
resulting gamma-ray flux, following Scott and Siverts- 
son. We then calculate the lower bound on the UMCH 
halo fraction which would result from detection of an 
UCMH by Fermi. We also calculate the upper bound 
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which would result if no UCMHs are detected. In Sec. lIIII 
we outline the calculation of the amplitude of the density 
contrast (c.f. Ref. Q), and translate the potential con- 
straints on the abundance of UCMHs into constraints on 
the power spectrum of the primordial curvature pertur- 
bation. We conclude with discussion in Sec. IIVI 



II. 



UCMH FORMATION AND GAMMA-RAY 
EMISSION 



Ricotti & Gould [6^ find that a density perturbation 
with amplitude at horizon crossing 5 > 10~^ will grow 
sufficiently during radiation domination that it collapses 
at z > 1000, seeding the formation a UCMH which then 
grows via spherical infall. It has been argued that PBHs 
can also seed the formation of minihalos [111, [l^ , and 
the resulting gamma-ray emission (assuming that the re- 
mainder of the dark matter is in the form of WIMPs) 
leads to constraints on the abundance of PBHs [l3[ , how- 
ever we do not pursue that possibility here. 

At matter-radiation equality the DM mass within a 
UCMH forming region, Af(zoq), is given by 01 



1 + Z, 
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where = iluu/^m — 0.834 [jyj is the dark matter 
fraction and M-a{zi) = {A7r/3)pH~^ is the horizon mass 
at redshift Zi corresponding to the epoch when the scale 
of interest entered the horizon. After matter-radiation 
equality the UCMH mass, Mucmh(z), grows, due to ra- 
dial infall of matter, as 
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Following Scott & Sivertsson [7| we assume that UCMHs 
stop growing at z w 10 as the onset of structure formation 
prevents further matter infall. Using the constancy of the 
entropy, s = (/^sO^r^, and the the radiation density, p = 
(7r^/30)(7*r'*, where (/^.s is the number of entropy degrees 
of freedom and gi, the number of relativistic degrees of 
freedom and T the temperature, the horizon mass can be 
written as 

/ eq \ 1/2 / rp \ 2 
MH(T)=MH(Teq)(^^j (^^j . (3) 

1 /3 

Using T oc g^^s (1 tlie horizon mass as a function 

of redshift is given by 



Mh(zO = MH(Zeq) ( ^) ^ (^i^^^ 
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where we have taken g^^s ~ g*- 

The UCMH dark matter density profile is given, in the 
radial infall model, by [1, 0| 



PvcMnir, z) 



3/xMucmh(^) 



167ri? 



UCMH 



(z)r^ 



(5) 



where RucMniz) is the radius of the UCMH at redshift 
z, given by 



^-Rucmh(z) 



pc 



= 0.019 



1000 

1 + z 



/ Mucmh(^:) 



(6) 



where Mq is the mass of the sun. 

Baryonic infall may lead to adiabatic contraction of 
the UCMH density profile [l3]. Scott & Sivertsson con- 
sidered a variable fraction of the total UCMH mass con- 
densing to form a constant density baryonic core. The 
dark matter density in the centre of the halo does not 
rise significantly and hence the change in the resulting 
gamma-ray flux is relatively small 0. Given the uncer- 
tainties in the calculation we therefore do not consider 
adiabatic contraction. 

WIMP annihilation reduces the density in the inner 
regions of the UCMH. We use the standard estimate of 
the maximum density, Pmax, 0j El] 



Pn 
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where is the WIMP mass, (av) the thermally aver- 
aged product of the WIMP annihilation cross-section and 
speed, to « 13.7 Gyr [ij the current age of the Universe 
and we take the UCMH formation time as ti{z = Zcq) « 
77 kyr [l^ . The UCMH present day density profile is thus 
given by pucMnir) = min {p^ax , Pucmh(»^, = 10)}, 
where pucmh(?', z — 10) is given by eq. 

The gamma-ray flux above a threshold energy _Eth, 
^-y{Eth), from WIMP annihflation within an UCMH at 
a distance d from the Earth can be written as 



^astro ^particle 

2d2 



The particle physics term, $particic7 is given by 

dNf 



particle 



/ 



dE 
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(9) 



where cr/ is the annihilation cross-section and dNf/dE 
the differential photon yield of the /th annihilation chan- 
nel. We use DarkSUSY [17|] to carry out a scan of the pa- 
rameter space of the Minimal Supersymmetric Standard 
Model and compute ^particle for the sets of parameters 
that are compatible with accelerator bounds and produce 
a present day DM density compatible with the WMAP 
measurement of the DM density. When calculating the 



^ This is true for dark matter in the form of standard WIMPs, with 
the canonical annihilation cross-section deduced from the mea- 
sured dark matter density. Motivated by recent electron data, 
Scott and Sivertsson also considered a model with enhanced an- 
nihilation cross-section. In that case WIMP annihilation leads 
to a larger constant density core and adiabatic contraction then 
has a larger effect. 
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lower limit on the halo fraction of UCMHs which would 
arise from a detection by Fermi we use the largest value of 
^particle obtained from the DarkSUSY scan. Conversely 
when calculating the upper limit which would result if 
no UCMHs are detected we use the smallest value. The 
astrophysical factor, ^astro, is given by 

Castro = / r^plcMBir,z^W)dr. (10) 

The Fermi point source sensitivity above 100 MeV 
is [11 

$T,(100MeV) = 6 X 10~^cm"2g-i ^ (^ll-^ 

For a given UCMH mass, Mucmh(2 = 0), we determine 
the distance d within which a UCMH of this mass would 
be detectable at threshold sensitivity by Fermi. We then 
calculate the fraction of the Milky Way in the form of 
UCMHs if there is a single UCMH within this distance. 
This is the smallest UCMH halo fraction which could 
be detected by Fermi. To do this we assume that the 
fraction of the DM in the form of UCMHs is independent 
of position so that the local and global UCMH fractions 
are identical 
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FIG. 1. Constraints on the UCMH halo fraction, /ucmh, 
as a function of present day UCMH mass, Mucmh(z = 0). 
The solid line shows the lower bound on the halo fraction 
which would result from the detection of gamma-rays from 
an UCMH by Fermi. The dotted line shows the upper limit 
on the halo fraction if gamma-rays from UCMHs are not de- 
tected, assuming DM is in the form of WIMPs. 
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where /3dm,mw('') is the density profile of the Milky Way 
halo, fT.ucMH,Mw(?') the number density of UCMHs and 
My)M,mw{< d) the mass of DM within a sphere of radius 
d centred on the Earth. We assume a NFW [l^ density 
profile for the Milky Way with parameters as found by 
Klypin et al. 

Fig. [T] shows the lower limit on the UCMH halo frac- 
tion, as a function of UCMH mass, which would result 
from the detection of a single UCMH by Fermi at thresh- 
old sensitivity. It also shows the upper limit on the 
UCMH halo fraction if Fermi does not detect gamma- 
rays from UCMHs, assuming that the DM is in the 
form of self-annihilating WIMPs. More massive UCMHs 
have a larger gamma-ray flux ($astro oc Mucmh(^ — 0) 
roughly) and hence can be detected at a larger distance 
{d cc A/ucmh(^ = 0)1/2). Mucmh(^ = 0) < 10^ Mq, 
d ^ lOkpc so that Mdm,mw(< d) increases more rapidly 
than A/ucmh(^ = 0) resulting in a decreasing limit 
on the halo fraction as Mucmh(2: = 0) is increased. 
For more massive UCMHs d becomes significantly larger 
than the scale radius of the Milky Way halo, hence 
AfDM,Mw(< d) cx In [Mucmh(-2 = 0)] resulting in a sub- 
sequent increase in the limit on the halo fraction for 
Mucmh(^ = 0) > IO^Mq. 



III. CONSTRAINTS ON THE PRIMORDIAL 
CURVATURE PERTURBATION 

To translate the limits on the UCMH halo fraction into 
constraints on the primordial curvature perturbation, we 
need to relate the present day UCMH halo fraction to 
the primordial density perturbation distribution. The 
present day UCMH density, fiucMH, is related to the 
UCMH halo fraction, /ucmh, by eq. P^ . Assuming that 
UCMHs are not destroyed by dynamical processes during 
structure formation, the present UCMH density is related 
to the fraction of the Universe at horizon entry which is 
overdense enough to later form UCMHs, /3ucmh, by 

f^UCMH — ^UM — N PvCM}i{Mii{zi)) (13) 

As UCMHs are far more compact and dense than typi- 
cal DM halos they will be far less susceptible to disrup- 
tion. Our lower bounds are conservative; if UCMHs are 
destroyed, the initial abundance of UCMH forming per- 
turbations, and hence the amplitude of the primordial 
perturbations, will be under-estimated. The upper limit 
from non-detection would, however, be weakened. 

If the smoothed density contrast, in the comoving 
gauge, ShoriR), at horizon crossing {R = {aH)~^), is in 
the range 10"^ < <5hor(-R) < 1/3, the DM in the region 
will eventually collapse to form an UCMH . The hori- 
zon mass Mii{zi) is related to the smoothing scale, R, 
by 1211] 
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MH(Zi) = AfH(^cq)(fccqi?)^ 



9\ 



(14) 



where fcoq — 0.07(f2m^^)Mpc~^ is the wavenumber and 
MYi{zaq) = 1.3 X 10''^(r2m/i^)~^ g the horizon mass at 
matter-radiation equahty and we use gl^ « 3 and g\ ~ 
100 0. 

The fraction of the Universe in regions dense enough 
to eventually form UCMHs is given by Press-Schechter 
theory [12], 

/•1/3 

/3uCMH(AfH(^i)) =2 / P(4or(i?))d(5hor(i?), (15) 

where, assuming that the initial perturbations are Gaus- 
sian, the probability distribution of the smoothed density 
contrast, P{5\^or{R)), is given by (e.g. Ref. [23j]) 
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where a{R) is the mass variance 



a\R)^ / w\kR)V5{k,t)— , 
Jo k 



and W{kR) = exp {—k'^R?/2) is the Fourier transform of 
the window function used to smooth the density contrast, 
which we take to be Gaussian. The relationship between 
the present UCMH density and the mass variance is then 
given by 
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The constraints on the present day UCMH density can 
therefore be translated into constraints on the mass vari- 
ance by simply inverting this expression. The final in- 
gredient required to calculate constraints on the power 
spectrum of the primordial curvature perturbation is the 
relationship between the density contrast in the comov- 
ing gauge and the primordial curvature perturbation, TZ. 
This has recently been calculated taking into account 
the evolution of perturbations prior to horizon entry (see 
Ref. 3 for details): 



6{k,t) 



4 / k 
71 \aH 



ji{k/V3aH)n, (19) 



where ji is a spherical Bessel function. The power spec- 
trum of the density contrast is therefore given by: 

MKt) = ^ fA-Vjf(k/V3aH)rTzik) . (20) 



Substituting this into eq. (|T7)) . and setting R = [aH) ^ , 
gives 

<AR) = Y X (kRfj!{kR/V3) 

dk 

xexp{-k^R^)Vn{k)^ . (21) 



This integral is dominated by scales k ^ ko = 1/R. 
Following Ref. :4] in the context of slow-roll inflation 
models we can assume that the power spectrum is con- 
stant over these scales, Vji^k) = V-niko)- Relaxing this 
assumption and using a power law power spectrum with 
spectral index in the range consistent with slow-roll in- 
flation, 0.9 < n{ko) < 1.1, leads to changes of order 3% 
in the power spectrum limits. Using eqs. ([T^ and (PT|) 
we can translate the UCMH abundance constraints into 
constraints on the amplitude of the spectrum of the cur- 
vature perturbation. For each UCMH mass considered 
we take the pivot point, ko, to correspond to the length 
scale of the perturbation (see eq. (HH)) which eventually 
forms the UCMH, ko = 1/R. 

Fig. [2] shows the constraints on the power spectrum 
of the primordial curvature perturbation. The poten- 
tial lower limit on the power spectrum which would arise 
from the detection of gamma-rays by Fermi from a sin- 



(17) gle UCMH is of the order V-r > 10" 
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10 
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scales k ~ 10^ — 10^ Mpc^^. If gamma-ray emission from 
UCMHs is not observed, an upper limit can be placed 
on the power spectrum of the primordial curvature per- 
turbation, of the order V-r ^ 10^^'^ — 10^^ on scales 
fc ~ 10^ - 10^Mpc~\ Constraints for lar ger wavenum- 
bers than those shown in Fig. [5] result in /ucmh ^ 1 and 
so are not considered. The lower bound based on a detec- 
tion at Fermi threshold sensitivity is a conservative limit 
(provided that the effects of adiabatic contraction are in- 
significant). The upper limit from non-detection relies 
on several assumptions, however, most significantly that 
the DM is in the form of WIMPs and that significant dis- 
ruption to UCMHs does not occur. If multiple UCMHs 
were detected by Fermi (or ACTs), or the fiux was sig- 
nificantly above the detection threshold, then this would 
imply a larger UCMH halo fraction, and hence the lower 
limits on the power spectrum of the primordial curvature 
perturbation would be stronger. 



IV. CONCLUSIONS 

Large amplitude (lO"'^ — 10^^) density perturbations 
can seed the formation of dense UCMHs Q . If DM is in 
the form of WIMPs, WIMP annihilation within UCMHs 
may produce a detectable fiux of gamma-rays [6, 7]. We 
have investigated the implications of detection, or non- 
detection, of gamma-ray emission from UCMHs by Fermi 
for the power spectrum of the primordial curvature per- 
turbation. We find that detection by Fermi at threshold 
sensitivity would place a lower limit on the UCMH halo 
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FIG. 2. Limits on tiie power spectrum of the primordial cur- 
vature perturbation as a function of comoving wavenumber 
(in units of Mpc~^). The sohd line shows the potential lower 
bound on the power spectrum resulting from the detection 
of gamma-rays from an UCMH by Fermi at threshold sensi- 
tivity. The dotted line shows the upper limit on the power 
spectrum obtained if gamma-rays from UCMHs are not de- 
tected by Fermi, assuming DM is in the form of WIMPs and 
UCMHs are not disrupted during structure formation. 



fraction. The smallest detectable UCMH halo fraction, 
/uCMH > 10-^ is for MucMH IO^Mq. 

If Fermi does not detect gamma-ray emission from 
UCMHs then, assuming DM is in the form of WIMPs, 
this would place an upper limit on the UCMH halo 
fraction. The limit is tightest, /ucmh 10^^, for 
MucMH ^ IO^Mq. The resulting potential upper limit 
on the power spectrum of the primordial curvature per- 
turbation, assuming UCMHs are not disrupted during 
structure formation, would be Vji ^ 10^^ ''' — 10^^ on 
scales k ^ 10^ — 10^ Mpc^"'^. This upper bound is signifi- 
cantly stronger than those from primordial black hole for- 
mation, Vtz < 10^^ — 10"^ 4], and would hence provide a 
tighter constraint on models of inflation (c.f. Ref. [1,1^). 
It does, however, rely on the assumptions that DM is 
in the form of WIMPs and UCMHs are not disrupted 
during the formation of the Milky Way halo. 
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